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Алгоритми обчислення спектрав сигналив 
в системах обробки зображень 


В статье рассматриваются алгоритмы вычисления спектров действительных сигналов, предназначенные 
для использования в системах цифровой обработки двумерных сигналов. Получен в удобном для 
практического использования виде алгоритм вычисления спектра действительной последовательности с 
помощью применения процедуры быстрого преобразования Фурье размером, равным половине длины 
реализации сигнала. Приведены структурные схемы двух алгоритмов. 
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Алгоритмы вычисления спектров сигналов в системах обработки изображений 


Введение 


Изображения, получаемые при картографировании поверхности Земли и дна 
Мирового океана радиолокационными и гидроакустическими методами, сейсмограм- 
мы, снимки, получаемые с помощью электронных микроскопов, представляют собой 
записи двумерных действительных сигналов. 

Для обработки таких сигналов используют двумерные системы цифровой обра- 
ботки изображений. Учитывая свойство разделимости двумерного преобразования Фурье, 
для определения спектров изображений применяют одномерные алгоритмы быстрого 
преобразования Фурье (БФ), которые дают максимальную эффективность при вы- 
числении дискретного преобразования Фурье (ДПФ) комплексных дискретных после- 
довательностей [1-6]. При использовании стандартных алгоритмов для вычисления 
спектров действительных сигналов их эффективность существенно снижается. 

Свойства симметрии спектров Х(А) действительных последовательностей х(иТ) 


Ве Х(К) = Ве Х(М-К), — Шш Х(Ю=-Шш Х(М-Ю, 
(| =|Х(М-Ю) аго Х (К) =-аго Х(М-Ю), 


позволяют с помощью выполнения одной процедуры БИФ размером № обработать 
одновременно две действительные последовательности по № отсчетов каждая или, 
используя процедуру БИФ размером №2, вычислить ДИФ действительной последо- 
вательности в М отсчетов, выполняя дополнительно несложные операции предобра- 
ботки и послеобработки. Такие совмещенные алгоритмы существенно уменьшают 
требуемое количество операций при цифровой обработке изображений. 

Целью данной работы является получение в виде, представляющем практиче- 
скую ценность, алгоритма вычисления ДПФ одной вещественной последовательности 
х(иТ) длиной М с помощью применения процедуры БИПФ размером №2, а также раз- 
работка структурных схем двух алгоритмов, ориентированных на обработку дейст- 
вительных сигналов. 


Описание совмещенных алгоритмов 


Рассмотрим совмещенный алгоритм вычисления ДИФ двух вещественных по- 
следовательностей х(иТ) и У(пТ) длиной по М отсчетов с помощью однократного ис- 
пользования процедуры БИФ размером №, заключающийся в следующем [7]. 

Из двух действительных последовательностей х(иТ) и у(иТ) формируется одна 
комплексная <(иТ) следующим образом: 


<(иТ) = х(иТ)+ У(иТ) ‚ где п=0, [...,М-1. 


Тогда ДПФ новой последовательности будет равно: 


М— 2” 
2(к)= &@ет)+ у@туе —м Е где К=0...М№-1. 


п=0 


Комплексно-сопряженное значение ДПФ для последовательности (ИТ) в точке 
(№-К) будет равно: 


7 (№ —К) = У [х(ит)- (ит) е 


п=0 


2 
Г п(М-К 
т ) 


У [жат)- рат) в : 
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Учитывая, что х(иТ) и у(иТ) действительные последовательности, можно записать: 
№1 в 
Ю-+И’(М-К)=2УжЖите ^ =2Х() 
п=0 


№МЬ—1 .2п 
Ю-И(м-Ю = ит" =277. 
п=0 
Откуда получаются аналитические выражения для ДПФ искомых действитель- 
ных последовательностей: 
2+2 (М-Ю 10-2 (и-Ю) 
О | 
2 йа 

На основе приведенных соотношений получена структурная схема совмещен- 
ного алгоритма вычисления ДПФ двух действительных последовательностей, изо- 
браженная на рис. |. 

Операция предобработки представляет собой простое формирование без до- 
полнительных вычислений комплексной последовательности <(пТ) из двух действи- 
тельных х(иТ) и у(пТ), а операция послеобработки - это вычисление с помощью 
приведенных выше формул искомых спектров Х (Ё) и У(К) по полученным значениям 
И (К), требующее дополнительно выполнения одного комплексного сложения на один 
спектральный отсчет каждой последовательности. Остальные операции послеобра- 
ботки - тривиальные. 


х(к) 


х(иТ) 
Операция 
У(иТ) предобработки 


Операция 
послеобработки У (к) 


п=0, 1,..., М-1 ЮО, Пиз 


Рисунок 1 — Структурная схема алгоритма вычисления ДПФ двух вещественных 
последовательностей 


Эффективность применения алгоритма заключается в уменьшении количества 
умножений по сравнению со стандартным БИФ в два раза. 

Для получения в удобном для практического использования виде совмещенного 
алгоритма вычисления ДПФ одной вещественной последовательности х(иТ) длиной 
М№ с помощью применения процедуры БИФ размером №2 представим любую крат- 
ную двум последовательность х(иТ) в виде двух последовательностей х/(пТ) и х2(иТ) 
длиной №/2 отсчетов каждая, сформированных из четных и нечетных отсчетов: 


х(пТ)= х(2иТ), гдеп=0,1,..., 5-1; 


х, ("пТ)= М(2п-+1Г], где п= 0, 1...., 5-1. 


Тогда №-точечное ДИФ исходной последовательности х(иТ) можно записать в виде: 


та 


Хх = Ут" +5 Мол+ те 


п=0 п=0 


2 
— 7 (Оп) 
Лу ) 
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Или, используя обозначения для четных х/(иТ) и нечетных х2(иТ) последователь- 
ностей: 


Ш ве о Е и. 
Х (Е) = Ух, (ит № м + Ух, (пе ЕЕ 
п=0 п=0 


хх (+ хде" 


где Х/(К), Х>(К) - №2-точечные ДПФ последовательностей х/(иТ) и х2(иТ). 

Последнее соотношение показывает связь ДПФ любой М№-точечной последователь- 
ности с ДПФ двух №/2-точечных последовательностей, составленных из четных и не- 
четных отсчетов исходной последовательности. 

Теперь положим, что исходная последовательность х(иТ) является действитель- 
ной. Тогда последовательности х/(иТ) и х2(иТ) также действительные. 

Сформируем новую комплексную последовательность длиной №2 отсчетов из 
двух действительных по следующему правилу: 


а(пГ) = х(иГ)+ } х.(иТ), где п=0, 1,2... 5. 


ДПФ последовательности <1(пТ) будет равно: 


м вы 
2(к)= Ут) + рыит) | ИР . где К =0..---1. 
п=0 


Тогда в соответствии с рассмотренным выше первым совмещенным алгорит- 
мом вычисления ДИФ двух действительных последовательностей с помощью одно- 
кратного использования процедуры БИФ имеем: 


х ()- РО, (И №. х- А-В И-Ю 


2 2) 


Откуда получим алгоритм вычисления ДПФ М-точечной исходной 
действительной последовательности х(иТ) с помощью БПФ размером №/2 в удобной 
для практического использования форме: 


кону" ет а - И 
] 


Де Л, МИГ, 

Значения Х (К) при К = №/2... М - 1 находят, используя свойства симметрии спектров 
действительных последовательностей. 

Структурная схема алгоритма изображена на рис. 2. 


х(иТ) Операция Операция 


предобработки послеобработки 


п=0,1,..., №/2 —1 К = О зы =- 


Рисунок 2 — Структурная схема алгоритма вычисления ДПФ одной вещественной 
последовательности с помощью БИФ размером №/2 
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В этом случае операция предобработки — это простое формирование №/2-то- 
чечной комплексной последовательности из четных и нечетных отсчетов исходной 
действительной последовательности, не требующее дополнительных вычислений. 

Операция послеобработки выполняется по последней полученной формуле и 
требует дополнительно выполнения одного нетривиального комплексного умножения и 
трех комплексных сложений на каждый спектральный отсчет. 

Эффективность алгоритма по количеству умножений можно оценить следующим 
образом: 


м 
о =. 
к.у. М М М 
— 10, —+М№М 108,—+4 
й 5 > 5 5 > 5 


2105 › № 


При М = 1024 эффективность равна 1,5. 


Выводы 


Полученный в удобном для практического использования виде совмещенный 
алгоритм вычисления ДИФ действительной последовательности х(иТ) длиной М№ с 
помощью применения процедуры БФ размером №2 и приведенные структурные 
схемы двух совмещенных алгоритмов ориентированы на вычисление двумерного 
спектра сигналов в системах обработки изображений. 
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Тре за4у \аз аппед ю оаш а!еогИБил$ оф геа| зедиепсе ОРЕТ сотршаНоп Гог 
ргасиса[ зе. 

А$ Ше гези, фе еНесйуе а1еогибт о РЕТ сотршайоп о геа| Л-- 1еп2 зедиепсе 
Бу теапз оЁ №/2-епощ ЕЕТ ргоседиге уаз 4еуеоре4, ап Фе Боск ФЧаетат$ оЁ о 
таре а1еогИтл$, пиепдеа Тог (\о-дитепз1юопа| $1епа1$ зресиит$ сотрщайоп шт 
зузетл$ ог Ч1еИа1 ппазе ргосеззште \аз$ оМатеад. 


Статья поступила в редакцию 05.06.2012. 
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